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Solid-state phase transformations and melting of high-purity crystalline boron have been in situ and ex situ 
studied at pressures to 20 GPa in the 1500-2500 K temperature range where diffusion processes become fast 
and lead to formation of thermodynamically stable phases. The equihbrium phase diagram of boron has been 
constructed based on thermodynamic analysis of experimental and literature data. The high-temperature part 
of the diagram contains p-T domains of thermodynamic stability of rhombohedral P-B 106 , orthorhombic 
y-&28> pseudo-cubic (tetragonal) t'-B 5 2, and liquid boron (L). The positions of two triple points have been 
experimentally estimated, i.e. P-t'-L at ~ 8.0 GPa and ~ 2490 K; and fi-y-t' at ~ 9.6 GPa and ~ 2230 K. 
Finally, the proposed phase diagram explains all thermodynamic aspects of boron allotropy and significantly 
improves our understanding of the fifth element. 
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he phase diagram of boron and thermodynamic stability of boron allotropes remain of fundamental interest 
in condensed matter physics and chemistry for a very long time - since the first characterization of pure 
element's allotropes in 1950s 13 . Later the existence of multiple boron modifications has been reported (see 
recent reviews 2,3 ) but some of them never have been reproduced, e.g. P'- and P"-B 4 , t-B 50 5 , "HP form" of Wentorf 6 , 
"HPHT form" of t-B 192 7 , etc. At present time only five allotropes are generally accepted: rhombohedral ot-B 12 
(oc-phase) 8 and P-B 10 6 (P-phase) 9 , orthorhombic y-B 28 (y-phase) 10 , tetragonal t-B 19 2 U and t-B 52 (t-phase in this 
paper) 5,12-14 . t-B 52 has been proved to exist only very recentiy and its crystal structure has not been unambiguously 
established so far. Two more phases have been predicted using ab initio structural evolution algorithm 15 , i.e. 
orthorhombic o-B 52 16 , closely related to t-B 52 , and metallic boron with a-Ga crystal structure 10 . The interest to the 
high-pressure behavior of boron has recentiy raised due to the discovery of boron superconductivity at high 
pressure 17,18 ; unusual partially ionic character of some B-B bondings in y-boron 10 ; high-pressure synthesis of novel 
boron-rich compounds 19 26 that are refractory and chemically stable 27 29 , superhard 3032 , and even could have 
metallic conductivity 33 ; unusual pressure-induced behavior of boron-containing icosahedral 34,35 and layered 36 39 
structures; and prediction of the nonmetal-metal phase transition in boron at a pressure above 89 GPa 10,40 . 

The stability of boron allotropes has been intensively investigated during past years using ab initio calculations. 
It has been predicted that at ambient conditions a-B 12 and P-B 10 6 have similar static energies, but disordered 
P-B 106 is more stable at ambient pressure, due to its lower zero-point vibrational energy 41 . At pressures above 
2 GPa, denser a-B 12 should be more stable 10 . The ab initio analysis of stability of boron structures 10 showed that at 
pressures above 20 GPa the a-phase loses its stability and another phase, orthorhombic y-B 28 (confirmed 
experimentally 10 ), becomes stable. Finally, above 89 GPa transition of semiconductive y-phase into metallic 
one should occur 10 . However, this pressure domain has not been explored experimentally at high temperatures 
(to overcome a kinetic barrier), and the latter allotrope remains to be discovered. Thus, numerous theoretical 
predictions require rigorous experimental studies. 

The first attempt to analyze the high-pressure phase equilibria in boron was made in 2007 42 based on ab initio 
calculations for ot-B 12 and P-B 10 6 phases and some ambient pressure experimental data. However, the result 
contradicts the experimental data on boron melting under pressure 43-44 , i.e. dT m /dp is overestimated by a factor of 
2 (see Fig. lc). Moreover, the reported a equilibrium temperature (Tp /Qt ) at ambient pressure is lower by 
~ 200 K as compared with the maximal temperature of P-to-oc recrystallization in the presence of Pt melt 45 . 

The first p-T phase diagram of boron was proposed only in 2009 by Oganov et al. 10 and contains 5 allotropes 
(four experimentally confirmed forms a, p, y, t-B 192 , and hypothetical metallic one) and liquid boron. This 
diagram combined more extended ab initio and experimental data on structural stability and phase relationships, 
but still some points remained unclear. Though overall correct, this phase diagram contained an uncertainty 
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Figure 1 (a) Sequence of synchrotron energy-dispersive X-ray diffraction patterns taken in situ in the course of stepwise heating of a-boron at 5.5 GPa 
(E m d = 78.07 eV A). Asterisk (*) indicates the position of the escape line of boron nitride (hBN, capsule material), (b) Experimental data on a-to-P 
transformation in boron. Open and solid circles represent the a- and P-phases observed during our in situ experiments. The up triangle corresponds to the 
beginning of the a-to-P solid-state transformation", while the down triangle is the onset of recrystallization of a-B 12 from P-B 106 in the presence of Pt 
melt 50 . The dashed line corresponds to the equilibrium between ot-B 12 and P-B 106 that has been calculated in Ref.42. (c) Experimental data on boron 
melting. Solid and open squares represent the crystalline and liquid boron observed during our in situ experiments. The crossed squares correspond to the 
beginning of melting. The open triangles show the literature data: (2350 K, ambient pressure) 60 , (2370 K, ambient pressure) 61 and (2480 K, 7.7 GPa) 43 . 
The dashed line represents the melting curve calculated in Ref.42. 



related to the stability field of the tetragonal boron phase, which at 
that time was thought to be "HPHT t-B 192 " 7 . A second point requir- 
ing elaboration is that the equilibrium line between a- and P-boron 
was an estimate, rather than direct measurement (which would be 
complicated by kinetics) or calculation (which would be complicated 
by disorder in P-boron). Recrystallization of tx-B 12 from P-B 10 6 was 
experimentally observed at much higher temperatures at both ambi- 
ent (from melts containing Pt 45 ) and high (during solid-state trans- 
formation 12 ) pressures. Finally, t-B 52 has been recently obtained at 
ambient 12 and high 1214 pressures, and even recovered as a single 



phase 1214 . Although this tetragonal allotrope has been interpreted 
as a metastable one as compared to mysterious "HPHT t-B 192 ", our 
recent results showed that the latter has a crystal structure related in 
many aspects to t-B 52 phase, rather than to t-B 192 46 . 

Very recently pseudo-cubic t'-B 52 of the t-B 52 structural family has 
been discovered 46 . It was recovered after experiments at 20 GPa and 
2500 K, the highest temperature reported so far for formation of a 
tetragonal phase. Contrary to common low-density t-B 52 phase(s) 
and related compounds, pseudo-cubic allotrope is quite dense, very 
close to y-B 28 . This phase seems to be a good candidate for a HPHT 
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allotrope, instead of strongly distorted "t-B 192 structure". At such 
high temperature the diffusion processes are quite intense even at 
such high pressure as 20 GPa. The observation of reversible trans- 
formation in boron at HPHT conditions would be a strong support to 
equilibrium phase diagram. 

In the present work we have studied the high-temperature part of 
boron phase diagram using in situ and recovery high-pressure 
experiments (to 20 GPa and 2500 K), as well as thermodynamic 
analysis. The established equilibria between boron allotropes and 
liquid phase are self-consistent from the point of view of classical 
thermodynamics and adequately describe the experimental data 
obtained by independent research groups. 

Results 

Phase transformation of <x-B 12 into P-B 10 6- Our in situ studies of 
a-boron in the 2-6 GPa pressure range using MAX80 multianvil 
system and energy-dispersive synchrotron X-ray powder diffrac- 
tion at HASYLAB-DESY 47 have shown that the phase transforma- 
tion of a-B 12 into P-B 10 6 occurs at temperatures above 1600 K 
(Fig. la & b). It passes noticeably higher the stability line predicted 
in Ref. 42 and shows strongly non-monotone temperature depen- 
dence in the 1600-1800 K range. At higher temperatures, 2400- 
2500 K, all lines of P-B 106 disappear due to the melting (Fig. la & c). 

Two opinions are know in the literature on the thermodynamical 
stability of a-B 12 , i.e. (i) at low and moderate pressures it has its own 
domain in thep-Tphase diagram 42 , (ii) below 8 GPa P-B 106 phase is 
stable allover the p-T region up to the melting temperature 43,48 . Since 
the a ±=*P equilibrium temperature obtained from ab initio calcula- 
tions 42,49 is lower by ~ 800 K than the onset temperature of the 
a-to-P transformation (Fig. lb), one cannot make a clear conclusion 
whether the crystallization of a-boron under these conditions is 
equilibrium or kinetically limited. 

Very recently an attempt to study the vicinity of the a~P~y triple 
point has been made by Parakhonskiy et al. 50 . Although the data on 
the growth of single-crystal a-B 12 in metal systems is interesting by 
itself, the authors' interpretation of the results in terms of equilib- 
rium thermodynamics is quite contradictory. From one side, at tem- 
peratures of 1400-1600 K the numerous metastable phases 
crystallize 12,51 , especially, on the timescale of Ref.50. Moreover, the 
extrapolation of the a 4=5 P equilibrium line reported in Ref.50 down 
to 1 MPa contradicts the crystallization of a-B 12 from Pt melt at 
higher (at least, 200 K higher than in Ref.50) temperatures 45 , i.e. in 
the domain of stability of P-B 10 6 phase. At such low temperatures and 
for such rigid structures the formation enthalpy of a-B 12 from P-B 106 
(Aif a /p) is more important for the estimation of thermodynamic 
stability of a phase than the fact of single crystal growth, which can 
occur outside the domain of thermodynamic stability, especially, in 
the case of covalent cage structures (e.g. Si clathrates 52 ). Although 
Parakhonskiy et al. 50 believe that their approach is similar to the 
classical study of the graphite ±5 diamond equilibrium, the close (in 
contrast to graphite/diamond) structural relationships between 
boron allotropes 4,10 can result in the metastable crystal growth of 
boron according to the Ostwald rule of stages. Moreover, the ability 
of boron to host transition metals (Pt, in particular) in the structure 53 
was completely neglected by Parakhonskiy et al. 50 , while it can influ- 
ence the mutual stability of various structures and the crystallization 
order (for example, metastable nitrogen-doped t-B 52 was observed in 
the B-BN system already at 5 GPa 20,22,29 ). 

To perform thermodynamic analysis of relative stabilities of a-B 12 
and P-B 10 6, we fitted experimental heat capacities to the Holzapfel 
equation 54,55 (Fig. 2a, Tab. l 28,56 58 ). The estimates for different con- 
tributions into total Gibbs energy AG a/ p (solid curves) are presented 
in Fig. 2b as dashed curves: (1) thermal due to the difference in heat 
capacities, and (2) configurational due to non-zero configurational 
entropy of P-B 10 6 1,49 - Taking into account only thermal part would 
lead to astonishing result: the a-B 12 phase should be stable either at 



high temperature or in the whole temperature range with very nar- 
row (500-750 K) stability domain of P-B 106 . This fact contradicts to 
all known experimental data on phase relationships between these 
two boron allotropes and cannot be overcome by variation of AiT^/p 
value. Only including the configurational entropy (the lower estim- 
ate based on structural data from Ref.59 is ~ 0.5 R Jmol _1 K -1 ) 49 
allows us to obtain reasonable AG a /p dependencies (Fig. 2b, solid 
curves). The impact of configurational entropy on thermodynamic 
properties is remarkable only in the case of P-B 10 6 (standard state): 
a-B 12 and y-B 28 have zero configurational entropy, while that 
of pseudo-cubic t'-B 52 is "hidden" in the corresponding ASx/r 
(quasi-constant at high temperatures) values describing the X ±=5t' 
equilibria. 

Now, for the a ±5 P "equilibrium" line reported in Ref.50 the Aif a/ p 
value should be about —4.5 kj/mol (Fig. 2b) but even this "lower 
estimate" for AHo/p seems to be too high as compared to both the 
results of ab initio calculations 10,49 and experimental value of 
AiT amorp hous/p 57 - Thus, the experimental estimation of AH a/ p (as well 
as S con fig. of P-B 106 ) seems to be the crucial point for establishing 
topology of the boron phase diagram in the vicinity of the tx— p— y 
triple point and at lower temperatures. So far, the a *=>fi and a ±=*y 
equilibrium lines reported in Ref.50 may be only considered as over- 
estimated upper boundaries of the thermodynamic stability domain 
of a-B 12 . Moreover, the most recent experimental 1 and theoretical 10 
studies clearly indicate that the a fc?P equilibrium line crosses the 
pressure axis, in contrast to Ref.50. 

Fig. 2c shows the tentative low-temperature part of the boron 
phase diagram. The 1 equilibrium line intersects the pressure 
axis, similar to previous simulations 1,10 . The slope (dp/dT) a _p = 
AS a _p/AV a _n is positive and has a non-zero value at 0 K due to the 
configurational enthropy of P-B 106 . Our lowest estimate for high- 
temperature S con fi g of P-B 106 is ~ 0.5 R that is consistent with 
high-temperature instability of a-B 12 . According to Monte-Carlo 
simulations (see Ref.l and references therein), the low-temperature 
Sconfig. should be ~ 5 times smaller than that at high temperatures. 
Thus, at low temperatures the a P equilibrium line should have the 
slope (dp/dTVp ~ (0.1 R J mok 1 K" ')/((). 153 cm 3 ) 1,8,59 . As for the 
(dp/dT) a _ T it tends to 0 when T — » 0; and at very low temperatures 
the a *=jy equilibrium line should be parallel to the temperature axis. 
The low-temperature part of the diagram thus constructed seems to 
be the most thermodynamicaUy consistent with all available experi- 
mental data and ab initio calculations, and is very close to that prev- 
iously reported by Oganov et al 10 . 

Melting curve of P-B 10 6- The experimental points on the P-B 106 
melting, both our experimental data and the previously reported 
values 43,44,60,61 are presented in Fig. 3. The disappearance of the 
diffraction lines of crystalline boron is accompanied by appearance 
of a broad halo typical for a liquid phase unambiguously points to the 
melting (Fig. la). Even slight decrease of the melt temperature leads 
to crystallization of P-B 10 6, indicative of equilibrium melting point 
and not metastable melting. At 5.45 GPa the melting temperature 
was found to be 2440 K. This is an intermediate value between the 
reported melting temperature of 2350-2370 K 60,61 at ambient 
pressure and 2480 K at 7.7 GPa 43 . The experiments at pressures 
below 2.5 GPa have always led to blowouts and have not allowed 
establishing the low-pressure part of the melting curve. 

According to Ref.62, Aif L/ p = 50.2 kj/mol at ambient pressure. 
For our calculations, the melting temperature of boron at ambient 
pressure has been chosen as 2360 K (the mean value of 2350 K 60 and 
2370 K 61 ). Since A V L/ p is the value that strongly influences the slope 
of the melting curve, we have used it as a fitted parameter. The zero- 
pressure atomic volume of P-boron has been calculated by extrapola- 
tion of the data on thermal expansion 63 . A good fit has been obtained 
with AV L/ p of 6.7% of the molar volume of crystalline boron at 
meting temperature which is close to the ~ 5-10% volume change 
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Figure 2 | (a) Heat capacity of cc-B 12 and P-B 106 at ambient pressure. The triangles show the experimental data 56-58 , while solid lines represent the data fit 
to the Holzapfel equation 54,55 , (b) Estimate for the formation enthalpy (H a ° K — Hp 0K ) of a-B 12 at ambient pressure (P-B 106 was considered as a standard 
state). The black dashed line shows the difference between Gibbs energies of a-B 12 and P-B 106 due to the non-zero configurational entropy of fj-phase, i.e. 
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0.5R)" 9 . The green dashed line represents the difference between thermal contributions (due to the difference in heat capacities) 



into the Gibbs energy of a-B 12 and P-B 106 . Blue solid line indicates the zero energy level chosen as Hp 0K . The highest and lowest temperatures of the 
synthesis of well-distinguishable a-Bi 2 crystals allow one to suggest that T^p should be between 933 K 5 °and 1500 K 51 (red rectangle area). Two red curves 
show the corresponding Gibbs energies of a-B 12 (G a — Gp) for two different values of (H a ° K — Hp 0K ); while these enthalpy values give the temperature 
limits (933 K 50 and 1500 K 51 ) for the a ±?P equilibrium (G a (T) - Gp(T) = 0). The mean value of (H a 0K - Hp" K ) is therefore about -4.5 kj/mol if 
crystallization of a = B 12 occurs at equilibrium conditions between these two temperatures, (c) Tentative low-temperature part of the boron phase 
diagram. The oc fc?P equilibrium line (dashed) crosses the pressure axis, similar to previous simulations 1,10 . The slope (d/>/dT) a _p = AS a _p/AV a _p is 
positive and has a non-zero value at 0 K due to the configurational enthropy of P-B 106 . The (dp/dT) a _ r ^ 0 whenT^> 0; and at low temperatures the a tsy 
equilibrium line (doted) is parallel to the temperature axis, while at higher temperature it joins the triple point defined by the intersection of the calculated 
a fc*p (dashed) and experimental p fesy (solid) equilibrium lines. 



reported in Refs.61, 64, 65 (Tab. 2) and is noticeably lower than the 
20% value taken for the calculation of the melting curve in Ref.42. 
Thus, the experimental melting curve of P-boron up to 8 GPa maybe 
described by the linear equation 



T L/p = 2360 +15.5 p 



(1) 



Table 1 | The 

capacity 5 ^ 58 

Parameters* 


parameters of the Holzapfel equation 5455 for heat 
of a-Bi2 and p-Bio6 

a-rhombohedral boron p-rhombohedral boron 


0, K 


374(49) 


970(10) 


Co 


2.38613 


0.06106 


c, 


-0.33881 


0.12374 


c 2 


1.6714 


0.52491 


A 


0.1 15 


0.125 



Equilibrium line between P-Bi 0 6 and y-B 2 8- Our experiments at 
pressures above 10 GPa have shown that samples quenched from 
1600-2000 K contain only pure y-B 2 s, while the samples quenched 
from 8 GPa contain only P-B 106 . These results well agree with 
experimental data of other groups 13,50 , and the P *=sy equilibrium 
line that fits both our and literature data can be described by equation 



T Y/P = 



-4534 + 707.8 p 



(2) 



Equilibrium line between J>-B 2 8 and t'-B 52 . Above 10 GPa the 
temperature seems to have more impact on the recovered boron 
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Figure 3 | Phase diagram of boron. The symbols show the experimental data. The solid lines represent the equilibria between different allotropes 
obtained by fitting the experimental points, as well as the unique melting curve of t'-B 52 thermodynamically consistent with other phase equilibria. Down 
and up triangles show solid and melted P-Bi 06 . Squares and circles correspond to y-B 2 s and t'-B 5 2, respectively, recovered after quenching down to normal 
conditions. Small open symbols show literature data 13 . Color guide for symbols and structures: blue - P-B 106 *, red - y-B 2 8 10 > wine - t'-B 5 2 4 ', black- liquid B 
(structural units of amorphous boron 75 ). 



Table 2 


Thermodynamic data on phase transformations in boron at pressures up to 25 GPa and temperatures between 1 500 and 3000 K 


Y/X 


AHy/x (kJ/mol), ASy/x (kJ/mol K) 
or AVy/x (cmVmol) 


Fitting option 


Experimental values and/or ab initio predictions 


L/P 


AH = 50 

AV = 15.5 AH/2360.0 = 0.32 
AS = AH/2360.0 = 2.1-1 0" 2 


Fixed to experimental value 
Constrained to experimental melting 
curve and experimental AH|yp 


AH me „ = 50.2 43 - 62 , 48. 93 70 

V me i, = 4.98 64 - 65 to5.22 61 ; Vp = 4.76" (V 3 ook = 4.57 59 ) 
AV= 0.22 to 0.46 


y/P 


AH = 2.5 

AV= - 707.8 AH/4534.0 = - 0.35 
AS= - AH/4534.0 = - 5.5- 10" 4 


Fixed to ab initio prediction 
Constrained to experimental 
transformation curve and AH y /p 


AH 0K = 
AV300K, 


2.5 10 

, MPa = 4.25 10 - 4.57 59 = - 0.32 


t'/P 


AH = 8.8 

AV= - 160.2 AH/3755.7 = - 0.38 
AS = AH/3755.7 = 2.3-1 0" 3 


Fitted 

Constrained to experimental AV300K, 
transformation curve and fitted AH t yp 


, M Pa = 4.27 46 - 4.57 59 = -0.30 


t'/y 


AH = AHf/p - AH r/p = 6.3 
AV= 5.1 AH/2176.0 = 0.01 
AS = AH/2176.0 = 2.9- 10~ 3 


Adjusted to fitted AH t yp and AH y /p 
Constrained to experimental 

transformation curve and 

adjusted AH t y. ( , 


AV300K, 


, MPa = 4.27 46 - 4.25 10 = 0.02 


L/t' 


AH = AH L/p - AH t7p = 41.2 
AV= AV L/P - AV typ = 0.7 
AS = ASl/p — AS t yp = 1.9-10 2 


Defined by AG,yp and AG|yp 
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allotrope as compared to the pressure. At 15-20 GPa, the samples 
quenched from 2500 K contained only mixture of y-B 28 and pseudo- 
cubic t'-B 52 46 , while below 2200 K t-B 52 (as well as a-B 12 ) could be 
observed as intermediate phase prior to crystallization of y-B 28 in 
the stability domain of the latter 1213 . At such high temperatures the 
kinetic factors do not play a decisive role anymore due to the intense 
diffusion. Thus, above 10 GPa pseudo-cubic t'-B 52 is stable at high 
temperatures, while y-B 28 - at low and moderate temperatures. 

All these results allowed us to define the domain of thermodyn- 
amic stability of y-B 28 in the phase diagram (Fig. 3). The triple point 
between (3-B 10 6, y-B 28 and t'-B 52 should be located at ~ 9.6 GPa and 
~ 2230 K, while the equilibrium line between y-B 28 and t'-B 52 is 
described by the equation 

T t/y = 2176 + 5.1 p (3) 

Only dense t'-B 52 phase can explain such a low pressure slope of the 
t' <=*y equilibrium line, contrary to other known tetragonal phases, 
t-B 52 and t-B 192 46 . 

Equilibrium line between P-B 10 6 and t'-B 52 . The position of the 
|3-y-t' triple point, experimental data on phase stability of P-B 10 6 
and the lowest pressure of t'-B 28 formation (in some experiments 
this phase has been recovered at ~ 7.7 GPa) allow one to draw the |3 
*=*t' equilibrium line (Fig. 3) which has negative pressure slope and 
follows the equation 

2V/p = 3756- 160.2 p (4) 

This equilibrium line intersects the (3-B 106 melting curve at 
~ 8.0 GPa and ~ 2490 K that is the triple point between P-B 106 , 
t'-B 52 and liquid. 

Discussion 

Among all boron allotropes mentioned above, the experimental ther- 
modynamic values are known only for a-B 12 and (3-B loe , at least at 
ambient pressure. That is why it seems reasonable to use |3-Bio6 as 
standard state with known thermodynamic potentials, while for L, t' 
and y phases (X or Y) the corresponding Gibbs energy should be 
corrected by a value of AGx/p defined by three parameters, i.e. AJJx/p. 
AV x/ p and ASx/p, independent of pressure and temperature in the 
first approximation. These parameters, if not known from experi- 
ment, can be obtained by fitting the experimental p-T data to theor- 
etical isopotential X *=jY equilibrium lines (i.e. AG X / Y (p,T) = 
AGx/p(p,r) - AGy/p^.T) = 0). 

Since the majority of experimentally observed equilibria between 
boron allotropes can be presented by straight lines, we did not com- 
plicated our analysis with exact formulae containing integrals. For 
example, p - A Vx/p and JAVx/r dp differ only by ~ 5% (particularly, 
due to the close compressibilities of all boron allotropes) 66 69 . Thus, 
the fitted parameters should be interpreted as average values over p-T 
domain of experimental data. 

In order to describe the high-temperature part of the boron phase 
diagram, one should know at least 9 parameters e.g. 3 expressions for 
AGx/p each containing 3 parameters, i.e. Aifx/p, AVx/p and AS X /p- 
These dependences follow the equation 

AGx/p = AH x/ p+pAV x/ p-TAS x /p. (5) 

For example, if one have such parameters for X = L, t' and y, two 
other Gibbs energies AG t -/ y and AGi/t' will be related to them by 
simple general relationship 

AG X / Y = AG X/ /p — AG Y /p (6) 

and will give us thermodynamically consistent equilibrium lines 
(AGx/p (p>T) = 0). In order to establish these 9 parameters, one needs 
at least 9 independent experimental values characterizing those equi- 
libria. Each experimental equilibrium line gives us two parameters: 
AH/ AS and A WAS according to the equation 



T eq = Ail/AS+p AV/AS (7) 

that is a good approximation at high temperatures. 

Four experimental equilibrium lines described above, i.e. L ±=>p, 
P «=5 y, t' P and t' ±^y give us 8 independent parameters. The ninth 
one we have chosen as experimental value of enthalpy of boron 
melting Aif L/ p (see Tab. 2) quite reproducible in independent experi- 
ments at ambient pressure 43,62,70 . 

The fitted and fixed values of the thermodynamic parameters are 
given in Tab. 2 (data from Refs.10, 43, 46, 49, 59, 61-65, 70). One can 
observe a reasonable agreement of the fitted parameters with esti- 
mates for A V L/ p and A V t »/p based on volume measurements at high 
temperatures, as well as with ab initio predicted Aif y/ p value at 0 K. 
All this additionally indicates self-consistency and uniqueness of 
thermodynamic equations used, as well as real physical meaning of 
the fitted parameters. 

The thermodynamic parameters of all mentioned above equilibria 
give the unique melting line of t'-B 52 , i.e. 

r L /f = 2186 + 37.4 p (8) 

The parameters that define corresponding Gibbs energy are given in 
Tab. 2. 

Figure 3 shows the phase diagram that represents all equilibria 
described above. At low pressures, two rhombohedral phases are 
stable, a-B 12 and P-B 10 6- These allotropes have related structures 4 , 
rhombohedrally distorted fee packing, produced by individual B 12 
icosahedra at low temperatures, and more complicated clusters at 
high temperatures. At pressures above 10 GPa, two other phases 
become stable, namely, y-B 28 and t'-B 52 . Both of them have struc- 
tures similar to distorted NaCl (orthorhombic and tetragonal distor- 
tions, respectively) 10,46 . It should also be noted that the common 
feature of two high-temperature phases (P-B 10 6 and t'-B 52 ) is a cer- 
tain degree of intrinsic structural disorder (partially occupied 
Wickoff positions), which gives to them non-zero configurational 
entropy already at 0 K and, subsequently, results in their stability 
at high temperatures. 

Finally, the equilibrium phase diagram of boron has been con- 
structed at pressures up to 20 GPa and temperatures up to 2500 K 
(Fig. 3). It has been experimentally proved that at least four boron 
phases, i.e. P-B 106 , y-B 28 , t'-B 52 and liquid boron, have p-T domains 
of thermodynamic stability. Two triple points have been established, 
i.e. the first between P-B 10 6, t'-B 52 and liquid at ~ 8.0 GPa and 
~ 2490 K; and the second between P-B 106 , y-B 28 and t'-B 52 at 
~ 9.6 GPa and ~ 2230 K. Thermodynamic analysis leads to the 
unique melting curve of t'-B 52 that is consistent with other experi- 
mentally established phase equilibria. 

Methods 

Multianvil experiments. Quenching experiments to 20 GPa were performed using 
6-8 type large-volume multianvil systems with octahedral pressure assemblies at 
Laboratoire des Sciences des Precedes et des Materiaux (LSPM-CNRS) and 
Bayerisches Geoinstitut (BGI). The experimental details are described elsewhere 71 ' 72 . 
Pressure and temperature have been either directly measured (thermocouples and 
p-T phase transitions of reference materials) or estimated from previously obtained 
p-Tcalibration curves; in all cases, the uncertainties were estimated to be about 1 GPa 
and 50 K, respectively. Samples were gradually compressed to the desired pressure at 
ambient temperature, and then the temperature was increased with a rate of about 
300 K/min up to the desired value. After heating for 5-10 min, the samples were 
quenched by switching off the power and then slowly decompressed. 

In situ X-ray diffraction. In situ experiments to 7 GPa were carried out using 
multianvil X-ray system MAX80 at beamline F2.1, DORIS III (HASYLAB-DESY). 
The experimental setup has been described elsewhere 47 . Energy-dispersive X-ray 
diffraction data were collected on a Canberra solid state Ge-detector with fixed Bragg 
angle 0 — 4.555(3) using a white beam collimated down to 100 X 100 urn 2 . 

The sample temperature up to 2200 K was measured by a W3%Re-W25%Re 
thermocouple. The correction for the pressure effect on the thermocouple emf was 
made using the data of Li et al 73 . Above 2200 K the power - temperature calibration 
curve was linearly extrapolated to the high-temperature region (up to 2600 K). 
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Pressures at different temperatures were found from the p-V-T equation of state of 
highly ordered (P 3 — 0.98 ± 0.02) graphite-like hexagonal boron nitride 74 . 

The samples were gradually compressed to the required pressure at ambient 
temperature and then diffraction patterns were collected at the stepwise (—50 K) 
temperature increase. With the storage ring operating at 4.44 GeV and 150 ± 50 mA, 
diffraction patterns were collected for 1 min in real time. After heating, the samples 
were quenched by switching off the power, and then the pressure was slowly released 
down to ambient. 

Ex situ X-ray diffraction. The recovered samples were studied by conventional 
powder X-ray diffraction using G3000 TEXT (Inel) diffractometer (Bragg- Br entano 
geometry) employing CuKal radiation. The goniometer was aligned with high purity 
silicon (a = 5.431066 A) and the standard sample of LaB 6 {a = 4.15695 A). 
Synchrotron X-ray powder diffraction measurements (k = 1.10347 A) have been 
performed at beamline 1711, MAX II (MAX-lab); Debye-Scherrer geometry with 
rotating quartz capillary was used. Unit cell parameters, the size of the blocks of 
coherent scattering and strains were derived from the LeBail profile analysis 
performed using the GSAS program. 
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